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ABSTRACT 

Taking advantage of the impressive sensitivity of Spitzer to detect massive galaxies at high redshift, we study the 
mid-infrared environments of powerful, high-redshift radio galaxies at 1.2 < z < 3. Galaxy cluster member candidates 
were isolated using a single Spitzer/IKAC mid-infrared color criterion, [3.6] — [4.5] > —0.1 (AB), in the fields of 48 
radio galaxies at 1.2 < z < 3. Using a counts-in-cell analysis, we identify a field as overdense when 15 or more red 
IRAC sources are found within 1' (i.e., 0.5 Mpc at 1.2 < z < 3) of the radio galaxy to the 5a flux density limits of our 
IRAC data {/a. 5 = 13.4/iJy). We find that radio galaxies lie preferentially in medium to dense regions, with 73% of 
the targeted fields denser than average. Our (shallow) 120s data permit the rediscovery of previously known clusters 
and protoclusters associated with radio galaxies as well as the discovery of new promising galaxy cluster candidates 
at z > 1.2. 

Subject headings: galaxies: active - galaxies: clusters: general - galaxies: high redshift - infrared 
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1. INTRODUCTION 

The classical technique of finding distant galaxy clus- 
ters from the extended X-ray emission associated with 
their intracluster medium is very effective at finding clus- 
ters out to z ~ 1, but rapidly becomes insensitive for 
cluster searches beyond z ~ 1.2 (for a detailed review 
of X-ray selected clusters, see Rosati et al. 2002). The 
challenge of using X-ray observations to find the most 
distant clusters is primarily due to the (1 -I- z)** fading 
of the X-ray surface brightness. While X-ray selection 
is very effective at finding massive structures at mod- 
erate redshifts and large samples are expected from the 
all-sk y eROSITA soft X-ray telescope after it launches in 
2013 (|Cappelluti et al.ll2011l ). X-ray selection is unlikely 
to find the most distant clusters and proto-clusters. 

Another proven method to identify distant clusters is 
by searching wide-field imaging surveys for the red se- 
quence of early-type galaxies. Such red sequence searches 
were initially introduc ed using solely optical data (e.g., 
IGladders fc Yeel[2000( ). and have identified hundreds of 
clusters out to z ~ 1. More recently, the Spitzer Adap- 
tation of the Red-sequence Cluster Survey (SpARCS; 
iWilson et al.l [2009:) has pushed this method to redder 
wavelengths and thus to higher redshifts. SpARCS has 
identified clusters out to z 1.3. However, a weak- 
ness of this approach is that it requires the presence of a 
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well-formed red sequence of early-type galaxies. Such a 
sequence might not yet exist as we approach the epoch 
of cluster formation. 

Other fiel d studies such as th e IRAC Shallow Sur- 
vev fISCS: |Stanford et all [20051 : IBrodwin et all [200fil : 
lEisenhardt et al.l 120081) also use d the Infrared Array 
Camera (IRAC: [Fazio et al.l l2004f) onboard Spitzer to ex- 
pand the sample of galaxy clusters known at z > 1.2. 
In particular, the ISCS identifies clusters on the basis 
of photometric redshifts and does not require a red se- 
quence. IRAC is an extremely sensitive tool for find- 
ing massive galaxies at high redshift since their 4.5^m 
flux densities remain nearly constant at 0.7 < z < 
2.5 due to a n egative and favorable /c-correction (e.g., 
IE senhardt et al.. 20081 

PapovichI (|2008[ ) used a simple IRAC color crite- 



rion (see Section 4) to highlight overdensities of high- 
redshift gala xies in the .gp ztzer Wide- Infrared Extragalac- 
tic f SWIRE: [Lonsdale et al.ll2003i) survey and, in doing 
so, isolated one of the highest redshift clusters known to 
date, CIG J0218-0510 at z = 1.62 in the AMM-LSS field 
of SWIRE. Spectrosco pic follow-up confirmed 9 members 
(jPapovich et al.l [Ml . This same c luster was indepen- 
dently discovered by 'Tanak a et al.l (|2010D using photo- 
metric redshifts. This second team found two concen- 
trations of galaxies at Zph ^ 1.6, one of them associated 
with extended X-ray emission. They spectroscopically 
confirmed six members. 

These various IRAC surveys have the strong advantage 
of providing uniformly selected galaxy cluster samples 
over wide areas. Such uniform cluster samples are ben- 
eficial for a range of studies, including probing the for- 
mation epoch of the early-type galaxy po pulations (e.g., 
iMancone et al.l[20Tol : iRettura et al.ll2011l| and statistical 
probes of cosmological parameters (e.g.. iVikhlinin et"all 
l200l[StCT^rern][M(Th . 

However, clusters are rare objects and finding larger 
samples of massive high-redshift galaxy clusters would 
require field surveys even wider than the several tens 
of square-degrees which is the current state-of-the-art. 
Work has been done using targeted cluster searches, 
i.e., focusing on regions of the sky suspected to host 
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overdensities of galaxies. For example, powerful, high- 
redshift radio galaxie^ (HzRGs hereafter) are among 
the most massive galaxies known up to very high red- 
shift (Seymour ct al. 2007) and, as such, are suspected 
to lie preferentially in overdense regions. Studies of 
HzRG environments have revealed excesses of extremel y 
red objects (EROs; iStern et al.l 12001 iBest et al.ll2Q0l . 
line emitters detected through narrow-ban d imaging 
(jPentericci et all 120001 : iVenemans et al.ll2007l and refer- 
ences therein) and submillimetre galaxies ('Steven s et al.l 
|T003; Dc Breuck ct al. 2004; Grevc et al. 2007) in the 
fields of HzRGs out to z ^ 5. A range of complementary 
studies have made use of optical through mid-infrared 
observations to isolate candidate cluster members asso- 
ciated with radio galaxies (Kodama et al. 2007; Galametz 
et al. 2009, 2010a; Hatch et al. 2010; Mayo et al. 2012 in 
press) . 

Recently, iFalder et al.l (j2010t ) extended this cluster se- 
lection technique to a large sample of high-redshift active 
galactic nuclei (AGN) at z ~ 1. Their sample, observed 
at 3.6/im with IRAC, included both radio-loud and radio- 
quiet AGN. They found excesses of 3.6/im sources within 
300 kpc of the AGN, as well as evidence for a positive 
correlation between source density and radio power. 

In this paper, we use mid-infrared 3.6 and 4.5/im ob- 
servations of radio galaxies at 1.2 < z < 3 to study their 
environments and identify their possible association with 
high-redshift galaxy clusters. Note that we use the desig- 
nation 'cluster' to refer to both cluster- and protocluster- 
like structures. The paper is organized as follows. The 
HzRG sample and catalog extraction are presented in §2 
and §3. In §4, we introduce our IRAC selection criterion 
for high-redshift cluster member candidates. In §5, we 
present the counts-in-cell analysis used to isolate over- 
densities associated with HzRGs. We also present the 
spatial distribution of IRAC-selected sources in our most 
promising high-redshift cluster candidates in §6. We dis- 
cuss our z > 3 fields in §7 and present our conclusions 
in §8. Throughout, we assume a ACDM cosmology with 
Ho = 70 km s-i Mpc^^, n„i = 0.3 and r^A = 0.7. Mag- 
nitudes and colors are expressed in the AB photometric 
system unless stated otherwise. 

2. RADIO GALAXY SAMPLE AND IRAC DATA 

Our primary HzRG sample comes from the Spitzer 
High-Redshift Radio Galaxy survey (SHzRG; Seymour 
et al. 2007, De Breuck et ah 2010). Most of the sources 
were observed with IRAC during Cycle 1 (PID 3329; 
PI Stern) with typical integration times of 120s in all 
four IRAC bands (see Table [2]). Some fields were imaged 
more deeply using guaranteed time observations (see Ta- 
ble EL^eteilson tlwIRAC reductions can b e foun d 
in ISevmour etHI (|2007D and IDe Breuck et al.l (|2010D . 
Mayo et al. 2012 (in press) report on the MIPS 24/im 
environments of 64 of these fields. 

The initial SHzRG sample was designed to homoge- 
nously cover the -Lsghz radio luminosity - redshift plane. 
More recent work on the SHzRG sample, including this 
paper, instead use J^snnM Uz, which is a rnore is otropic 
measure of AGN power (|De Breuck et al.ll2010f ). This 

^ Following Seymour et al. (2007), we define a HzRG as a radio 
galaxy above a redshift of one with a restframe 3 GHz luminosity 
greater than 10^^ W Hz-^ 



leads to a slightly less uniform distribution in the 
-^sooMHz radio luminosity - redshift plane (see Fig. 2 of 
De Breuck et al. 2010). We will keep this in mind by 
using appropriate statistical tests later in the paper (see 
Section 5.3). 

In 2009, we obtained deeper 1600s IRAC observations 
of a sample of 10 radio galaxies (see Table[21 PID 60112; 
PI Hatch) at 3.6 and 4.5/im. Six of these fields were al- 
ready part of the initial shallower sample. In this paper, 
we make use of the deeper data w hen availab l e. Th ese 
deeper data were reduced as iuiSe vmour et al.l (j2007D . 

The goal of our study is to isolate galaxy clusters 
at high redshift using only mid-infrared data. As de- 
scribed below, the IRAC criterion applied in this work 
(see Section 4) is optimal for isolating galaxy structures 
at 1.2 < z < 3. We therefore concentrate our analysis 
on the 48 HzRG fields in this redshift range. We also 
analyze the 25 HzRG fields at z < 1.2 and 3 < z < 5.2 
from our sample in an identical manner and use these 
fields as a control sample. 

3. CATALOG EXTRACTION AND FLUX LIMITS 

In order to study the mid-infrared environments of our 
HzRG sample, we analyzed the available Spitzer data 
using standard techniques. 

We restricted the analysis to areas covered for at least 
60s in both IRAC channel 1 (3.6/im) and 2 (4.5/im). 
The central part (~ 1' x 1') exposure time is 120s. We 
checked the images by eye and flagged zones affected 
by artifacts (e.g., scattered light). Source extract ion 
was done using SExtractor (jBertin fc Arnoutill996D in 
dual image mode, using the 4.5/im frame as the de- 
tection image. W e used SExtractor parameters from 
iLacv et all ()2005D which have been optimized for anal- 
ysis of IRAC data, and measured photometry in 3" di- 
ameter apertures. These magnitudes were then corrected 
to total magnitudes using corrections determined by the 
IRAC instrument team (M. Lacy, private communica- 
tion). Specifically, the aperture corrections applied were 
1.68 and 1.81 for channels 1 and 2 respectively. 

Limiting flux densities for each image were determined 
from randomly placed 3" diameter apertures. For the 
SHzRG fields, the limiting flux densities were determined 
on the area covered by 120s exposure which extends be- 
yond a radius of 1' from the HzRG i.e., the size of the cell 
on which this work is conducted. To allow for uniform 
analysis of our large data set, we adopt a conservative 
flux density cut corresponding to the 5a depth of our 
shallowest data, i.e., /se = 11.0/.tJy ([3.6] = 21.3) and 
/4.5 = 13.4/i Jy ([4.5] = 21.m 

Recently, iMancone et al.l ()2010f ) derived the 3.6 and 
4.5/im luminosity functions of galaxy clusters out to 
z - 1.5. They found [3.6]* ~ 20.0 and [4.5]* ^ 20.1 
for clusters at z ~ 1.2. Similarlv. iStrazzullo et al.l ()2006l ) 
derived the -R's-ba nd luminosity function of clusters at 
z ~ 1.2. Using the lBruzual fc Char lot! ()2003D model of a 
2 Gyr-old single burst galaxy with an exponentially de- 
clining star formation history with r = 0.1, we find that 
[3.6] ^ 0.3 at z ~ 1.2 with a nearly c onstan t K— [3.6] 
color at 1 < z < 2. IStrazzullo et all (|2006l ) find that 
K; - 20.5 at z ~ 1.2, corresponding to [3.6]* 20.2, in 

Magnitudes in IRAC channels 1 and 2 are indicated by [3.6] 
and [4.5], respectively, and mAB = 23.9 — 2.5 log(/^/l/iJy). 
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agreement with the work of iMancone et ahl (|2010f) . To 
our adopted hmiting magnitude of [3.6] = 21.3, our shal- 
low 120s observations are therefore sensitive to galaxies 
on e magnitude fa i nter th an L* at z ^ 1.2. As shown 
in IMancone et al.l (|201ClD . m* should remain relatively 
constant out to z ~ 2 assuming cluster galaxies form at 
z > 2.5. 

4. COLOR SELECTION OF CLUSTER MEMBER 
CANDIDATES 

One of the main spectral features in the spectral en- 
ergy distributions (SEDs) of galaxies is the 1.6/im bump, 
caused by a minimum in the opacity of the H ~ ion 
which is present in the atmospheres of cool stars ()Johnl 
[19881) . This bump, seen in the SEDs of all normal 
galaxies, has often been u sed as an efficient photomet- 
ric redshift indicator (e.g., ISimpson fc Eisenhardtl[r999l : 
ISorba fc Sawick i 20101 and references therein). Fig. [T] 
shows the position of the 1.6/im bump for galaxies at 
z = 1 and z = 2 relative to the IRAC 3.6 and 4.5/im 
bands. The bump enters the IRAC bands at z ^ 1 and 
shifts beyond 3.6/im at z 1.2. This causes lower red- 
shift galaxies to have blue [3.6] — [4.5] colors, while galax- 
ies above z ~ 1.2 become red across these passbands. 

Combining a wide variety of composite stellar popu- 
lation models and spectroscopy of high-redshi ft galaxies 
in the GOODS-S and AEGIS DEEP2 fields, iPapovichI 
(pool confirmed that [3.6]-[4.5] > -0.1 (i.e., h.e/U.b < 
1.1) is very efficient at isolating z > 1.2 galaxies. They 
showed that only a minority population of0.2<z<0.5 
strongly star-forming galaxies with a strong warm dust 
contribution have such red colors and that about 80% of 
galaxies with [3.6] — [4.5] > —0.1 are at z > 1. Specif- 
ically^ based on rcdshifts from the DEEP2 redshift sur- 
vey, iPapo vich (2008) show that ^ 50% of galaxies at 
z > 1.1 have red IRAC colors, with this percentage in- 
creasing to ^ 75% for z > 1.2 galaxies and reaching 
~ 90% for z > 1.3 galaxies. In the later analysis we use 
this color cut to study the fields of HzRGs at z > 1.2, 
though we caution the reader than the efficiency of the 
selection criterion is diminished in the 1.2 < z < 1.3 
range. At any rate, only one HzRG from our sample, 
3C266 at z = 1.275, is in this redshift range. Finally, 
we note that, contrary to other high-redshift galaxy se- 
lection t echniques such as the BzK selection for z > 1.4 
galaxies (jPaddi et al.ll2000l ) or the near-infrared selection 
techniques for z > 1 .6 galaxies (jKajisawa et al. 2006; 
iGalametz et al.l[2010bl ). this single IRAC color criterion 
does not permit the segregation between galaxy types 
since it is based on a spectro-ph otometric prop erty of 
essentially all galaxy populations ()Papovichll2008[ ) . 

This red IRAC color criterion will identify a few 
additional astronomical populations, but none are ex- 
pected to be significant contaminations. Mos t stars have 
[3.6] — [4.5] ^ —0.5 (i.e., Vega color of zero). iStern et all 
(|2007f ) shows that only brown dwarfs cooler than spectral 
type T3 have such red colors ([3.6] - [4.5] > -0.07) and 
such objects will be rare in a fiux-limited survey at our 
depth. Extremely dusty stars, such as asymptotic giant 
branch (AGB) stars, w ill also have red IRAC colors (e.g., 
lEisenhardt et al.]|2010D . However, they are likewise ex- 
pected to be quite rare in a flux-limited survey primarily 
targeted at high Galactic latitudes. 

At all redshifts, most powerful AGN (> 95%) have 
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Fig. 1. — Spectral energy distribution of a 2 Gyr-old galaxy red- 
shifted to z = \ (dashed line) and 2 = 2 (solid line) assuming a 
rapid (r = 0.1 Gyr) exponentially declining star- formation history. 
The strongest feature in this wavelength range is the restframe 
1.6/im stellar bump. Transmission curves for IRAC channel 1 and 
2 are overplotted for reference. 

[3.6] - [4. 5] > -0.1 and wil l be selected by this IRAC 
criterion. IStern et al.l (l2005h measure a surface density 
of 275 AGN per deg^ in the Bootes field down to flux 
density limits of /s.g — 12.3^Jy and /4.5 = 15.4/.tJy. Us- 
ing IRAC catalogues available for the Bootes field (see 
Section 5.2 for details), we derive a density of ^ 28, 700 
red IRAC-selected sources per deg^ to the same depth. 
We therefore expect contamination by AGN to represent 
less than 1% of sources with [3.6] — [4.5] > —0.1 at the 
depth of our survey. 

According to models for passively evolving stellar pop- 
ulations formed at high-redshift, negative fc-corrections 
provide a nearly constant 4.5 fim flux density at z > 0.7 
and up to high redshifts. Fig.[2]shows the evolution of the 
4.5/im flux density and the [3.6] — [4.5] color as a function 
of redshift for various formation redshifts. We make use 
of the publicly av ailable model calculator EZ Ga0. We 
use the lMarastonI (l200a) and the Ch ariot fc Bruzual 2007 
— an update of iBruzual &: Chariot ()2003l ) — models for 
single stellar populations assuming solar metallicity and 
a Salpeter initial mass function. The 4.5/im flux density 
is normalized to match the observ ed ml of galaxy clus- 
ters at z - 0.7, [4.5]vcga = 16.75. IMancone etldl (|20100 
studied a large sample of clusters up to z ~ 1.5 to derive 
the 4.5/im luminosity function and found that the lumi- 
nosity function is consistent with a formation redshift 
z/ ^ 2 — 3. Fig.[2]shows that for z/ = 3, the 4.5/im flux 
density is relatively constant at z > 1. It even decreases 
(brightens) with redshift as one approaches the galaxy 
f ormation epoch. 

IMancone et al.l (|2010f ) models could not simultane- 

^ www. mancone. net/ezgal/model 
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ously match the low and high-redshift data; they found 
that m* becomes fainter at the highest redshifts probed, 
interpreted as a possible evidence for g alaxy assembly 
at z > 1.5. Note, however, that the IMancone et al.l 
uncertainties increase with redshift, with few of 
the highest redshift clusters spectroscopically confirmed, 
and larger systematic uncertainties in subtracting the 
foreground/background galaxy populations. Numerous 
authors have investigated the epoch of early-type galaxy 
assembly, and the results span a wide range of forma- 
tion redshifts (e.g., van Dokkum 2008; Eisenhardt et a]] 
I2010f ). A larger sample of high-redshift galaxy clusters is 
clearly needed to further investigate this question. 

In the following, we assume a uniform [3.6] — [4.5] > 
—0.1 color criterion to select high-redshift candidates in 
the fields of 48 radio galaxies at 1.2 < z < 3. We consider 
sources detected to the 5cr depth of our SHzRG sample 
([4.5] = 21.1). We include sources with [3.6] - [4.5] > 
—0.1 with a 3.6/im magnitude fainter than the limits of 
our SHzRG survey. For such sources, we uniformly assign 
a lower limit of [3.6] = 21.3. Sources selected by these 
criteria are referred to as 'red IRAC-selected sources' or 
simply as 'IRAC-selected sources' in the rest of the paper. 

As seen in Fig. [21 our selection criteria should robustly 
identify clusters across this whole redshift range, assum- 
ing cluster galaxy formation and assembly is at similar 
or higher redshift. The aim of the current work is to test 
this hypothesis. Current literature only studies evolved 
populations in clusters out to z ^ 1.3, with less certain 
results for cluster candidates at higher redshifts. Exten- 
sive literature from multiple observing programs infer the 
cluster formation epoch is at higher redshift, with results 
ranging over the 2 < zj^ < 5 range. Given the current 
state of knowledge, the high redshift end of our analysis 
is somewhat arbitrary; we could have been more conser- 
vative and chosen a lower redshift cut-off, or, in principle, 
we could have considered a higher redshift cut-off. If we 
do not find rich fields around our highest redshift HzRGs, 
multiple explanations are possible. It could mean that 
clusters do not exist at those redshifts, it could be due 
to small number statistics, or it could imply issues with 
the simple evolutionary model we are testing in Fig. [21 
a model that is remarkably effective at describing clus- 
ter data at z < 1.3. In particular, galaxy mergers might 
cause an evolution of m* with redshift such that we be- 
come less sensitive to clusters at the high-redshift end of 
our sample. Such a result would be of interest. 



Fig. 2. — Evolution of 4.5^im flux density (top panel) and [3.6] — 
[4.5] color (bottom) versus redshift for single stellar populations 
with a range of formation redshift, zj = 2 (blue), 3 (green) and 5 
(red). Rapid (r = 0.1 Gyr) exponentially declining star-formation 
history models, from Maraston et al. 2005 (M05; solid line) and 
Chariot & Bruzual 2007 (CB07; dashed), were generated with EZ 
Gal (assuming [4.5]vGga = 16.75 at 2 ~ 0.7; see text for details). 
At z > 1, [4.5] is relatively independent of redshift out to the 
formation epoch. The dotted lines indicate the [4.5] depth (top) 
and the [3.6] — [4.5] > —0.1 color criterion (bottom) adopted in this 
work. 



5. COUNTS-IN-CELL ANALYSIS 
5.1. Counts-in-cell in SWIRE 

We use a counts-in-cell analysis to identify overden- 
sities of IRAC red galaxies associated with HzRGs. 
As a reference field, we analyze the SWIRE (SWIRE; 
iLonsdale et al.l 120031 ) survey. With its typical 120s ex- 
posure, the SWIRE survey reaches slightly deeper than 
our HzRG data, primarily because it exclusively targets 
low-background regions of the sky. SWIRE covers ^ 50 
deg^, split into six independent fields which mitigates the 
effects of cosmic variance. 

To ensure consistency with our analysis of the radio 
galaxy fields, we build our own catalogues for all six 
SWIRE fields. As described in Section 3, we use the 
4.5/im image for source detection and derive photometry 
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Fig. 3. — Distribution of the surface density of [3.6] — [4.5] > —0.1 sources in 20,000 cells of 1' radius randomly generated in the six 
SWIRE fields to the flux density limits of our sample (dotted histogram). The corresponding number of sources per cell is indicated on 
the right axis. A Gaussian fit of the lower part of the distribution is shown by the red curve and the 2a deviation above the mean of this 
distribution by the horizontal dotted line. Each radio galaxy is placed at the vertical axis location corresponding to the surface density of 
IRAC-selected sources found within 1' from the radio galaxy. The solid histogram shows the percentage of radio galaxy fields for a given 
density. Two examples of known high-redshift galaxy clusters (CLG J0218-0510 and CI J1449-I-0856) are also indicated by the gray labels 
(see §5.2). 

density regions) by a Gaussian (iteratively clipping 
at 2a). The best fit values give a Gaussian centered 
at (N) — 2.80 galaxies per arcmin^ 8.8 sources 
per cell) with crjv = 0.97. We define an overdensity 
of IRAC-selected sources when a cell is denser than 
N + 2(7 jv ^ 4.74 galaxies per arcmin^ — i.e., when at 
least 15 sources are found within a cell of 1' radius. This 
criterion, for which visual inspection typically registers 
interesting spatial segregation, is matched by only 9% 
of the cells in SWIRE. 



from 3" diame ter apertures, cor rected to total magni- 
tudes following ILacv et all (|2005D . 

At the depth of our survey, we detect the bright end 
of the galaxy luminosity function, which preferentially 
populates cluster cores. We therefore adopt a cell size of 
1' radius, corresponding to ~ 0.5 JVIpc at 1 < z < 3. This 
distance matches typical cluster r200 size s for z > 1.3 
mid-i nfrared selected galaxy clusters (e.g.. iWilson et all 
l2009l) . 

We generate 20,000 independent (i.e., non- 
overlapping) circular cells of 1' radius, randomly 
distributed over the six SWIRE fields. Fig. |3| shows the 
histogram of the surface density of red IRAC-selected 
sources per cell. The horizontal axis shows the percent- 
age of cells corresponding to a specific surface density. 
The shape of the distribution is a Gaussian strongly 
skewed towards higher density cells. We fit the lower 
half of the distribution (i.e., the distribution of the lower 



5.2. Tests on known z > 1 galaxy clusters 

Before applying our newly defined criterion, we test 
its efficiency on published clusters at z > 1 ob- 
served with IRAC. We first study two of the highest 
redshift c lusters known to date: CIG J0218-051 at 
z = 1.62 fPap ovich et al.l |2010^ 'Tanaka et al.'l2010[) and 
CI J1449-t-0856 at z = 2.07 (fGobat et al..,20l£ 
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CIG J0218-0510 is located in the XMM field of SWIRE. 
We select sources with [3.6] — [4.5] > —0.1 in the field of 
CIG J0218-0510 using the catalogues we derived from the 
SWIRE data (see Section 5). CI J1 449+0856 was di scov- 
ered in the CNOC 1447+09 field (lYee et al.l[2000!) and 
was also first identified as an overdensity of red IRAC- 
selected galaxies. To date, 1 1 cluster members h ave been 
spectroscopically confirmed (jCobat et al.ll20lil ). We re- 
trieved the IRAC 3.6 and 4.5^m reduced post-BCD im- 
ages of the CNOC 1447-1-09 field from the Spitzer archive 
(P.I. Fazio) and derived our own catalogue in an identical 
manner as for the HzRG targets. 

For both cluster fields, we select sources with [3.6] — 
[4.5] > —0.1. Fig. |4] shows the spatial distribution of 
these red sources in the 6' x 6' field around CIG J0218- 
0510 (top) and CI J1449-f 0856 (bottom panel). Large red 
symbols account for sources detected within the limits 
applied to our SHzRG sample and smaller gray dots ac- 
count for sources within the 5a limits of SWIRE (/a g = 
3.7fi3y and /4.5 = 5.4^Jy) and CNOC 1447-h09 (/a.e = 
2.3^Jy and /4.5 = 3.8^Jy). We also show the spec- 
troscopically confirmed cluster members (black squares) 
for C IG J0218-0510 (Paoovich et al. 2010; Tanaka et S\ 
|2010[ ) and CI J1449+0856 (R. Gobat, private communi- 
cation). All of the confirmed members detected at 3.6 
and 4.5/um have [3.6] — [4.5] < —0.1 though some are 
below the 5a limits plotted in Fig. |4l 

Nineteen red IRAC-selected galaxies are found within 
1' of the assigned center of CIG J0218-0510 and 15 for 
CI J1449-I-0856. Both clusters are clearly seen as com- 
pact concentrations of red galaxies and are selected as 
overdense fields by our IRAC criteria (see also Fig. [31 
gray labels). 

An overdensity of EROs was also recently detected in 
the fi eld of the quasar 3C270.1 at z = 1.53 (iHaas et al.l 
|2009| ) suggesting that the quasar is associated with 
a high-redshift galaxy structure. We retrieved the 
IRAC post-BCD 3.6 and 4.5/.tm images of this field 
from the Spitzer archive and extracted the correspond- 
ing 4.5/im selected catalogue. Twenty-four red IRAC- 
selected sources are found within 1' of the quasar. This 
strengthens past results that the quasar is part of a high- 
redshift clus ter. 

Recently, lEisenhardt et al.l ()2008l) identified a large 
sample of 335 mid-infrared selected galaxy clusters in 
the Bootes field, including five clusters spectroscopically 
confirmed at 1.2 < z < 1.5. The Bootes field was cov- 
ered by the Spitze r Deep, Wide-Field Survey (SDWFS; 
lAshbv et all [20091 ) ■ an IRAC survey of 8.5 deg^ with 5a 
depths of 22.9 and 22.4 in 3.6 and 4.5/xm, respectively. 
We make use of the publicly available 4.5/im-selected 
catalogue of the SDWFS Data Release 1 to derive the 
number of IRAC-selected sources within 1' of the cluster 
centers. All five z > 1.2 confirmed galaxy clusters have 
15 or more IRAC-selected sources in the studied cell and 
would be selected as overdense fields by our IRAC cri- 
teria. Of the 40 spectroscopically confirmed members of 
these clusters at z > 1.2 detected at 4.5/im (Stanford et 
al. 2005, Brodwin et al. 2006, Eisenhardt et al. 2008), 
25 (63%) have [3.6] - [4.5] > -0.1. Considering only 
clusters members at z > 1.3, all 21 galaxies have colors 
consistent with the red IRAC criterion within the photo- 
metric errors. Considering the z < 1.2 clusters reported 
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Fig. 4. — Spatial distribution of sources with [3.6] — [4.5] > -0.1 
in the fields of two spectroscopically confirmed galaxy clusters: 
CIG J0218.3-0510 (z = 1.62; top) and CI J1449+0856 (z = 2.07; 
bottom; the solid line marks the limit of the area covered by 
Spitzer). Small gray circles indicate sources down to the 5a depth 
of the original IRAC data while larger red symbols indicate sources 
detected to the flux limits of our SHzRG data. Dotted circles illus- 
trate 0.5, 1 and 2 Mpc radii at 2 = 1.62 and z = 2.07. A 1' radius 
cell (~ 0.5Mpc) centered on the adopted cluster center is shown 
by the red dashed circle. Spectroscopically conflrmed members are 
indicated by black squares. North is up; East is to the left. 

in Eisenhardt et al. (2008), only the cluster at z = 1.161 
would be selected as an overdense field by our red IRAC 
selection criterion. These results confirm that our IRAC 
criterion begins to become effective at identifying galaxy 
clusters at z > 1.2 and is quite robust at identifying rich 
structures at z > 1.3. 

A large number of high-redshift galaxy cluster candi- 
dates were also recently found in SWIRE by the SpARCS 
survey. We further test our selection criterion on spectro- 
scopically confirmed SpARCS clusters and recover their 
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two highest redshift galaxy clu sters (to date) : SpARCS 
J003550-431224 [z = 1.34; Wil son et al .112003) — a very 
dense structure with 35 IRAC-selected sources found 
within 1' of the assigned cl uster center — and S pARCS 
J161037+552417 (z = 1.21: lDemarco et al.llMol) with 19 
red IRAC-selected sources. We note that the two z < 1.2 
clusters published in Demarco et al. 2010 (at z = 0.871 
and z = 1.161) were not recovered. 

5.3. Radio galaxy fields 

We now determine the number of red IRAC-selected 
sources within 1' of each radio galaxy. Each field is shown 
in Fig. |3] at the position corresponding to the number of 
IRAC-selected sources found within 1' of the radio galax- 
ies. The density for each field is also reported in Tables 2 
and 3 (last column) . The solid histogram shows the per- 
centage of radio galaxy fields corresponding to a given 
density. We find that radio galaxies preferentially lie 
in medium to dense regions, with 73% (±12%) of the 
targeted fields denser than the mean of the Gaussian 
fit to the SWIRE fields. Eleven fields i.e. 23% (±7%) 
are found with more than 15 red IRAC-selected sources 
within 1' of the central radio galaxy. We conduct a two- 
sided Kolmogorov-Smirnov statistic between the distri- 
butions of densities of IRAC-selected sources in SWIRE 
and in the HzRG cells and find a probability of only 1.3% 
that the two distributions are drawn from the same dis- 
tribution. 

From the 20, 000 cells we produced in SWIRE, we gen- 
erated 100, 000 random samples of 48 cells and derived 
the fraction of these sub-samples that have at least 11 
2(T-overdense cells. We find that only 0.3% of the sub- 
samples have a fraction of overdense cells comparable to 
our radio galaxy sample, confirming that radio galaxies 
lie preferentially in dense regions. 

We discuss our 11 overdense fields in more detail in 
Section 6. At our bright magnitude limits, we expect to 
isolate high-redshift galaxy structures that contain many 
bright (M* -|- 1 or brighter) galaxies in their cores. We do 
not rule out the possibility that the other radio galaxies 
may be associated with less massive and/or less compact 
galaxy structures. 

The sample of radio galaxies covers a wide range of 
redshift as well as a wide range of radio luminosity. We 
next consider a possible dependance of overdensity with 
these two parameters. Fig.[S]and Fig.lHlshow the counts- 
in-cell around the radio galaxies versus redshift and radio 
luminosity, respectively. We adopt Poissoni an errors for 
the counts-in-cell and use the lGehrelsl (fT986l) smaU num- 
ber approximation for Poisson distribution. The inset 
on top of each figure shows the percentage of overdense 
fields compared to the total number of fields per bin of 
redshift or radio luminosity. 

No significant correlation is observed between over- 
density and redshift (Fig. [5]). As far as a correlation 
with radio power is concerned (Fig. [5]), we note that 35% 
(±13%) of radio galaxies with Lsoomhz > lO^^'^ W Hz"! 
are found in overdense fields, compared to 23% of the 
full sample. This could indicate a positive correlation 
between overden sity and radio power , a result suggested 
by p ast studies ()Milev fc De Breuckl [20081: iFalder et all 
I2010D . However, the Spearman rank correlation coeffi- 
cient between overdensity and radio power is as low as 
0.071. As mentioned in Section 2, there may be a small 
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Fig. 5. — Number of [3.6] — [4.5] > -0.1 sources within 1' of the 
radio galaxies (and within the 5cr 4.5^m Umit of our sample) as a 
function of redshift. Overdense environments are marked in red. 
The horizontal red line shows the median value (8.8) for a 1' radius 
cell in SWIRE. The top panel shows the percentage of overdense 
fields relative to the total number of fields per bin of redshift. 

degeneracy between L500MHZ and redshift in our sample. 
A study of the Spearman partial rank correlation coef- 
ficient — which takes the co rrelation of tw o variables 
with a third related variable (|]VIacklinlll982t ) — still do 
not show any significant correlation between overdensity 
and radio luminosity. This result is consistent with the 
absence of correlation observed between overdensity and 
the density of 24/im JVIIPS-selected sources in the same 
HzRG sample (IVIayo et al. 2012 in press). We conclude 
that a larger sample, reaching radio luminosities of sev- 
eral orders of magnitude lower, is required to study (and 
quantify) a possible correlation. 

5.4. Analysis of HzRG fields with deeper IRAC data 

We also examine fields with at least 1600s IRAC ex- 
posures. We adopt correspondingly deeper cuts for the 
deeper sample (74.5 = 3.4/iJy; [4.5] = 22.6). Similar to 
the shallower analysis, we select sources using the crite- 
rion [3.6] — [4.5] > —0.1. For sources below the deeper 
3.6/im fiux density limits, we uniformly assign a lower 
limit [3.6] = 23.0. The adopted depth corresponds to 
^ M* + 3. We derive densities of red IRAC-selected 
sources within 1' of the radio galaxies to these deeper 
fiux limits (Table [T]). We find results consistent with 
our analysis at shallower cuts, recovering PKS1138-262 
(z = 2.156), JVIRC0156-252 (z = 2.016) and MRC0350- 
279 (z = 1.900) as the densest fields; compared to the 
average surface density of red IRAC-selected sources in 
the other radio galaxy fields, the environments of these 
three radio galaxies are at least 1.5 times richer. 

6. CONFIRMED CLUSTERS AND NEW CLUSTER 
CANDIDATES 

Eleven radio galaxy fields present significant overden- 
sities of red sources. Fig.[7]shows the spatial distribution 
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Fig. 6. — Number of [3.6] - [4.5] > -0.1 sources within 1' of the 
radio galaxies as a function of radio luminosity. Symbols are as 
in Fig. [5] We note that the majority of our cluster candidates are 
found around HzRGs with LsooMHz > 10^^'^ W Hz~^ suggesting 
that more powerful radio galaxies reside in denser environments. 



TABLE 1 

Surface density of IRAC-selected sources for the deeper 

FIELDS. 



HzRG 




Density 






(arcmin~^) 


PKS1138-262 (2 = 


= 2.156) 


18.1 ± 2.4 


MRC0156-252 {z 


= 2.016) 


17.8 ±2.4 


MRC0350-279 (z 


= 1.900) 


17.5 ± 2.4 


MRC1017-220 (z 


= 1.768) 


14.3 ±2.1 


MRC2104-242 (z 


= 2.491) 


14.0 ±2.1 


MRC2139-292 (z 


= 2.550) 


13.7 ±2.1 


MG2308+0336 (2 


= 2.457) 


12.1 ± 2.0 


USS1425-148 (2 = 


= 2.349) 


12.1 ± 2.0 


MRC0324-228 {2 


= 1.894) 


11.5 ± 1.9 


LBDS53W002 (2 


= 2.393) 


10.8 ± 1.9 


MRC0406-244 (2 


= 2.427) 


10.8 ± 1.9 


MRC1324-262 (2 


= 2.280) 


8.6 ± 1.7 



of the IRAC-selected sources for these 11 cluster candi- 
dates. Small gray and large red dots indicate red IRAC- 
selected sources detected within the 3cr limits of each in- 
dividual field and to the depth of our survey, respectively. 
These fields present a very diverse spatial distribution of 
red sources. Some fields show clumps (e.g., USS1425- 
148, 1VIRC0350-279) while others show more filamentary 
structures (e.g., MRC0114-211, 3C470, PKS1138-262). 

The fields of JVIRC0114-211 {z = 1.41), 3C470 (z = 
1.65), JVIRC2224-273 (z = 1.68), MRC0350-279 [z = 
1.90), MRC2139-292 (z = 2.55) and MRC2025-218 (z = 
2.63) have never been studied for overdensities in the 
past. They show several clumps of red sources whose 
association with the targeted radio galaxy still need to 
be confirmed with spectroscopy. We next briefly discuss 
each of these cluster candidates in order of increasing 
redshift. 



6.1. MRCOII4-2II fz = 1.4i; 

The environment of this radio galaxy had never been 
studied in the past. This field is the lowest redshift but 
also the strongest galaxy cluster candidate in the HzRG 
sample, overdense at a 6(j level compared to SWIRE. 
Only 0.03% of the cells in our counts-in-cells analysis are 
found with such density. The distribution of the IRAC- 
selected red sources shows a clear spatial segregation. 
The bright red galaxies lie preferentially in two N/S elon- 
gated structures with the radio galaxy embedded in the 
southern one. The northern clump is 2' (~ 1 JVIpc) north 
of the radio galaxy. Additional groups of red galaxies are 
also observed within 1 JVIpc of the radio position, suggest- 
ing that JVIRCOl 14-211 is part of a complex large-scale 
structure. Furthermore, five red IRAC-selected sources 
(including the radio galaxy itself) are found in the im- 
mediate vicinity (within 15" ^ 125 kpc) of the HzRG, 
corresponding to a density about 9 times higher than av- 
erage. No other radio galaxy in our sample shows such 
a dense environment. 

6.2. 7C1756+6520 (z ^ IA2) 

iGalametz et al.l ((2009) identified this HzRG as resid- 
ing in a candidate high-redshift cluster based on BzK 
photometry. This cluster was recent ly confirmed in 
s pect roscopy with 21 cluster members fcalametz et al.l 
2Q10b). The spatial distribution of the IRAC-selected 
sources is clumpy, with a main compact clump of bright 
galaxies 1' SE of the radio galaxy. This location coincides 
with t he spectroscopic a lly con firmed sub-group of galax- 
ies in IGalametz et all (|2010bf ) at z = 1.437. Seven of 
the spectroscopically confirmed members were observed 
and detected bot h at 3.6 and 4.5/x m. Only one (ID 
CI 1756.20 in Galametz et al.l [20101)0 has an unexpect- 
edly blue IRAC color ([3.6] - [4.5] = -0.4). The other 
six galaxies have [3.6] — [4.5] > —0.1. 

6.3. 3C470 (z = 1.65; and MRC2224-273 (z = IM) 

These fields present very inhomogenous distributions 
of IRAC-selected sources. In both cases, the radio galax- 
ies are located (in projection) at the edges of the galaxy 
overdensities. Spectroscopic follow-up of the structures 
is required to confirm the association of high-redshift 
galaxy structures with the radio galaxies. 

6.4. MRC0350-279 (z = 1.90; 

This is the third densest field in our sample. Two 
very compact clumps of red sources arc found within 1' 
and aligned with JVIRC0350-279 in the direction NW/SE, 
forming an intriguing symmetrical structure. 

6.5. MRC0156-252 (z = 2.02; 

IGalametz et al.l ()2010aD presented a study of the sur- 
roundings of 1VERC0156-252. In that work, we used a 
purely near-infrared color criterion to select z > 2 galax- 
ies in the field and found that the radio galaxy lies in 
an overdense region of both early-type and star-forming 
candidates. The majority of the IRAC-selected sources 
within 1' of 1VIRC0156-252 are found north of the radio 
galaxy, a result also found for the ne ar-infrared selected 
cluster candidates in iGalametz et al.1 (,2010a,). 
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6.6. PKS1138-262 (z = 2.16) 

This radio galaxy (also referred to as the "Spider 
Web" galaxy) is associated with a well studied high- 
redshift galaxy protocluster. Several intense programs 
have studied this field, identifying a di versity of cluster 
members: Lya emitters, Ha emitters (jPentericci et al.l 
dpOO; Kurk ct al. 2004a), or candidate cluster members: 
EROs (Kurk ct al. 2004b), Distant Red Galaxies (DRGs; 
[Kpdama ct al. 2007). To date, only two massive red 
sources have been spectrosco picaUy confirrned to be as- 
sociated with PKS1138-262 (jPohertv et al.ll2010[ ): both 
are red IRAC-selected sources with [3.6] — [4.5] > 0.1. 
Similar to DRGs selected in this field, the (brightest) 
IRAC-selected galaxies are preferentially found along a 
3'-long E / W filamentary structure (corresponding to the 
radio axis) in which the radio galaxy is embedded. 

6.7. USSI425-I48 (z = 2.35; 

The environment of this radio galaxy was recently 
studie d using near-infrared observations. iHatch et al.l 
(|2011| ) selected high-redshift galaxies using near-infrared 
color criteria in the field of six radio galaxies at z ^ 2.4. 
USS1425-148 was the most overdense field of their sample 
and suspected to contain a high-redshift galaxy cluster 
at z > 2. Fig. [7] suggests that the radio galaxy has sev- 
eral nearby companions. A second clump of sources is 
found ~ 1' north of USS1425-148. 

We note that MRC1324-262 (z = 2 28) and MIIC0406 - 
244 {z = 2.427) were also part of the IHatch et al.l ()2011l ) 
sample and were not found to lie in particularly overdense 
environments, a result consistent with the present mid- 
infrared analysis of these fields. 

6.8. 4C23.56 (z = 2.48^ 

4C23.56 is also known to lie in a n overdense region of 
galaxies. More than a decade ago, iKnopp fc ChambersI 
(|1997f l foimd an overdensity of red galaxies (J — if > 4, 
Vega) in the immediate surroundi ngs of 4C23. 56. A com - 
plementary analysis was made bv lKaiisawa et all ()2006f ) 
who found that 4C23.56 was the densest field out of their 
sample of six radio galaxies at 2.3 < z < 2.6 in terms 
of near-infrared cluster galaxy candidates. An excess of 
candidate Ha emitters has also re cently been detected 
in the field bv iTanaka et al.l ()201lD . which also reported 
on the first spectroscopicaUy confirmed companions for 
4C23.56, three Ha emitters at z = 2.49. 

6.9. MRC2139-292 (z = 2.55^ and MRC2025-218 
(z = 2.63; 

These fields are our two highest redshift candidates. 
The field of MRC2139-292 presents a very clear and com- 
pact concentration of red galaxies centered on the radio 
galaxy. Though MRC2025-218 is just above our thresh- 
hold for identifying rich IRAC environments, we note 
that Mayo et al. 2012 (in press) find that it resides in a 
very rich field in the MIPS 24/im bandpass, overdense at 

7(7 level compared to control fields. 

7. THE MID-INFRARED ENVIRONMENTS OF Z > 3 
RADIO GALAXIES 

The initial sample of HzRGs also comprised 21 galaxies 
at z > 3 that were analyzed in an identical manner to the 
main sample. Some of these radio galaxies are known to 



be within spectroscopicaUy confirmed high-redshift pro- 
toclusters (e.g. USS0943-242, TNJ1338-1942, MRC0316- 
257, TNJ0924-2201; Venemans et al. 2007) though none 
were recovered by the present mid-infrared analysis. This 
suggests evolution in the luminosity function of high red- 
shift galaxy clusters. As seen in Fig. 2, the simplest sin- 
gle stellar population models predict relatively constant 
[4.5] flux densities at z > 1 for a wide range of formation 
redshifts, as well as a brightening of the galaxies as we 
approach the formation epoch. The fact that we do not 
recover these known rich (proto-)clusters at z > 3 implies 
that this simple model is incorrect, most likely because 
the massive ellipticals that comprise the clusters are still 
in process of forming in a hierarchical manner. Deeper 
IRAC data from a large sample of clusters and cluster 
candidates out to z ~ 3 should test this hypothesis. We 
are, in fact, amassing such a data set with the Clus- 
ters Around Radio-Loud AGN (CARLA) Warm Spitzer 
snaphot program (Galametz et al. in prep.). 

The radio galaxy 4C41.17 (z = 3.79) lies in an over- 
dense region of IRAC-selected sources. A compact clump 
of IRAC-selected sources is observed, slightly offset 1' 
south of the radio galaxy. Ivison ct al. (2000) also re- 
ported an overdensity of sub-millimeter galaxies in this 
field, though spectroscopic follow-up has f ailed to con- 
firm their association with the radio galaxy (jCreve et al.l 
l2007f) . We suspect that there might be a foreground 
galaxy structure in this field. 

8. CONCLUSIONS 

We identified overdensities in the fields of powerful 
HzRGs at z > 1.2. We studied the distribution of sources 
with [3.6] — [4.5] > —0.1 using a counts- in-cell analysis, 
defining a region to be overdense when denser (at the 2a 
level) than the SWIRE average, i.e., when 15 or more red 
IRAC-selected galaxies (with 74.5 > 13.4/iJy) are found 
in a cell of 1' radius. We identified 11 radio galaxies with 
at least 15 IRAC-selected candidates within 1'. These 11 
radio galaxy fields, representing 23% of our sample of 48 
HzRGs, are the most promising galaxy cluster candidates 
from our analysis. The large fraction of overdense fields 
clearly indicates that radio galaxies, on average, reside in 
rich enviroments. Randomly selecting similarly analyzed 
samples of 48 control fields from the SWIRE survey, we 
find such a high fraction of dense fields in only 0.3% of 
the samples. 

Five of these fields have been studied in the past and 
are known or suspected to contain a (proto-) cluster at 
high-redshift. Three have been spectroscopicaUy con- 
firmed to be associated with the radio galaxy. This 
clearly demonstrates the strength of IRAC to select high- 
redshift galaxy structures. Our relatively shallow data 
have permitted the recovery of known clusters in very 
short integration times (120s); past ground-based stud- 
ies have required hours of optical/near- infrared observa- 
tions. 

We tested for correlations of environment with radio 
galaxy properties and found no significant trends with 
either redshift or radio luminosity. We note however 
that a detailed analysis of evolution of the cluster pop- 
ulation with redshift will require an even larger sample 
of high-redshift clusters as well as a better understand- 
ing of the luminosity function for clusters at z > 1. No 
statistically significant trend with radio luminosity was 
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Fig. 7. — Spatial distribution of [3.6] — [4.5] > —0.1 galaxies around the radio galaxies with a 2a excess of IRAC-selected sources compared 
to SWIRE (i.e., with more than 15 IRAC-selected sources within 1'). We plot sources down to the Sit limits of each image (small gray 
dots) and to the uniform (5(t) limits of our full sample (larger red dots). Contours of the combined IRACl+2 studied zone are illustrated 
by the dashed lines. Concentric dotted circles account for distances of 0.5, 1 and 2 Mpc at the redshift of the targeted radio galaxy (yellow 
star). The red dashed circle shows the cell of 1' radius centered on the HzRG and used in the count-in-cell analysis in Section 5. North is 
up, East is to the left. 



observed. Nevertheless, the majority of our cluster can- 
didates are found in the fields of more luminous radio 
galaxies (Lsoomhz > lO^^ "^ W Hz"!). 

This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy under a contract with NASA. We are very grateful to 
Mark Brodwin and Peter Eisenhardt for having provided 
information on the Bootes cluster sample mentioned in 
this paper and to Conor Mancone for providing his useful 
EZ Gal Model Generator and valuable help on models. 
We also thank the anonymous referee for his/her very 
useful comments. 
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TABLE 2 

The shallow HzRG sample (in redshift order). 





R A. 


Dec 


2^ 


log(L500MHz)^ 












(W Hz"-"-) 


( T/^vn 1 Tl ^ ^ 
( 0,1 Lllilli ) 


3C356.0 


17:24:19.0 


50:57:40.30 


1.079 


28.35 


2.5 ± 0.9 


MRC0037— 258 


00:39:56.4 


—25:34:31.01 


1.100 


27.72 


3.5 ± 1.1 


3C368.0 


18:05:06.3 


11:01:33.00 


1.132 


28.52 


3.2 ± 1.0 


6C0058+495 


01:01:18.9 


49:50:12.29 


1.173 


27.33 


3.8 ± 1.1 


3C266 


11:45:43.4 


49:46:08.24 


1.275 


28.54 


3.8 ± 1.1 


MRC021 1—256 


02:13:30.5 


—25:25:21.00 


1.300 


27.78 


2.5 ± 0.9 


MRC0114-211 


01:16:51.4 


-20:52:06.71 


1.410 


28.66 


8.9 ± 1.7 


7C1756+6520 


17:57:05.4 


65:19:53.11 


1.4156 


27.40 


6.4 ± 1.4 


7C1751+6809 


17:50:49.9 


68:08:25.93 


1.540 


27.46 


2.5 ± 0.9 


3C470 


23:58:35.3 


44:04:38.87 


1.653 


28.79 


5.7 ± 1.4 


MRC2224-273 


22:27:43.3 


-27:05:01.71 


1.679 


27.52 


4.8 ± 1.2 


6C0132+330 


01:35:30.4 


33:17:00.82 


1.710 


27.64 


4.1 ± 1.1 


3C239 


10:11:45.4 


46:28:19.75 


1.781 


29.00 


2.2 ± 0.8 


3C294.0 


14:06:44.0 


34:11:25.00 


1.786 


28.96 


2.9 ± 1.0 


7C1805+6332 


18:05:56.9 


63:33:13.14 


1.840 


27.78 


4.5 ± 1.2 


6CE0820+3642 


08:23:48.1 


36:32:46.42 


1.860 


28.28 


1.6 ± 0.7 


6CE0905+3955 


09:08:16.9 


39:43:26.00 


1.883 


28.17 


3.5 ± 1.0 


6CE0901+3551 


09:04:32.4 


35:39:03.23 


1.910 


28.19 


3.2 ± 1.0 


MRC0152-209 


01:54:55.8 


—20:40:26.30 


1.920 


28.20 


3.8 ± 1.1 


MRC2048— 272 


20:51:03.6 


—27:03:02.53 


2.060 


28.72 


4.5 ± 1.2 


5C7.269 


08:28:38.8 


25:28:27.10 


2.218 


27.82 


2.2 ± 0.8 


4C40.36 


18:10:55.7 


40:45:24.01 


2.265 


28.79 


2.9 ± 1.0 


TXS0211— 122 


02:14:17.4 


— 11:58:46.00 


2.340 


28.48 


4.1 ± 1.1 


USS1707+105 


17:10:06.5 


10:31:06.00 


2.349 


28.63 


1.6 ± 0.7 


USS1410— 001 


14:13:15.1 


—00:22:59.70 


2.363 


28.41 


1.3 ± 0.6 


6C0930+389 


09:33:06.9 


38:41:50.14 


2.395 


28.41 


2.2 ± 0.8 


3C257 


11:23:09.2 


05:30:19.47 


2.474 


29.16 


3.2 ± 1.0 


4C23.56 


21:07:14.8 


23:31:45.00 


2.483 


28.93 


6.4 ± 1.4 


USS1558— 003 


16:01:17.3 


—00:28:48.00 


2.527 


28.82 


3.8 ± 1.1 


WNJ1115+5016 


11:15:06.9 


50:16:23.92 


2.540 


27.82 


1.3 ± 0.6 


Ubb0828+193 


08:30:53.4 


"1 o 1 n "1 r\r\ 

19:13:16.00 


2.572 


28.44 


3.2 ± 1.0 


PKS0529— 549 


05:30:25.2 


—54:54:22.00 


2.575 


29.16 


2.2 ± 0.8 


MRC2025— 218 


20:27:59.5 


—21:40:56.90 


2.630 


28.74 


5.7 ± 1.3 


U&D2202+128 


22:05:14.1 


13:05:33.50 


2.706 


28.54 


2.2 ± 0.8 


MG1019+0534 


10:19:33.4 


05:34:34.80 


2.765 


28.57 


3.2 ± 1.0 


4C24.28 


13:48:14.8 


24:15:52.00 


2.879 


29.05 


2.9 ± 1.0 


4C28.58 


23:51:59.2 


or* 1 r\ r\c\ 

29:10:28.99 


2.891 


28.91 


3.8 ± 1.1 


Usb0943— 242 


09:45:32.7 


— 24:28:49.65 


2.923 


28.62 


2.9 ± 1.0 


WNJ0747+3654 


07:47:29.4 


36:54:38.09 


2.992 


28.14 


2.5 ± 0.9 


B3J2330+3927 


23:30:24.9 


39:27:12.02 


3.086 


28.33 


1.3 ± 0.6 


WNJ0617+5012 


06:17:39.4 


50:12:55.40 


3.153 


28.02 


3.9 ± 1.1 


MRC0251— 273 


02:53:16.7 


—27:09:13.03 


3.160 


28.54 


2.5 ± 0.9 


WNJ1123+3141 


11:23:55.9 


31:41:26.14 


3.217 


28.51 


2.9 ± 1.0 


6C1232+39 


12:35:04.8 


39:25:38.91 


3.220 


28.93 


2.5 ± 0.9 


TNJ0205+2242 


02:05:10.7 


22:42:50.40 


3.506 


28.46 


1.6 ± 0.7 


i IN JU12i + lo2U 


Ui:21:4z. / 


13:20:58.00 


3.516 


28.49 


3.2 it 1.0 


TXJ1908+7220 


19:08:23.7 


72:20:11.82 


3.53 


29.12 


5.4 ± 1.3 


USS1243+036 


12:45:38.4 


03:23:20.70 


3.570 


29.23 


4.1 ± 1.1 


WNJ1911+6342 


19:11:49.6 


63:42:09.60 


3.590 


28.14 


1.9 ±0.8 


MG2144+1928 


21:44:07.5 


19:29:14.60 


3.592 


29.08 


3.5 ± 1.1 


6C0032+412 


00:34:53.1 


41:31:31.50 


3.670 


28.75 


1.9 ±0.8 


4C60.07 


05:12:54.8 


60:30:52.01 


3.788 


29.20 


4.1 ± 1.1 


TNJ2007-1316 


20:07:53.3 


-13:16:43.62 


3.840 


29.13 


3.2 ± 1.0 


8C1435+635 


14:36:37.1 


63:19:14.00 


4.250 


29.40 


2.5 ±0.9 


TNJ0924-2201 


09:24:19.9 


-22:01:41.00 


5.195 


29.51 


1.0 ±0.6 



Our analysis focusses on fields around HzRGs with 1.2 < z < 3. The other fields are used as control fields. 
Radio luminosities from lDe Breuck et al.l 1 12010 ). 
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TABLE 3 

The deeper Spitzer sample (in redshift order). 



Name 


R A 


UCC. 


z ^ 


log(-i'500MHz)'' 


irv,/\v-/ ijjxp. lime 




Density 




J ZUUU j 






(\\J TT-y— 1\ 

vv nz ) 


\p) 




(arcrQin ^) 


LBDS53W091 


17:22:32.7 


50:06:01.94 


1.552 


27.04 


900 


65 


3.8 ± 1.1 


MRC1017-220 


10:19:49.0 


-22:19:58.03 


1.768 


27.94 


1600 


60112 


4.5 ± 1.2 


MRC0324— 228 


03:27:04.4 


—22:39:42.60 


1.894 


28.49 


1600 


60112 


2.9 ± 1.0 


MHCU35U— 279 


03:52:51.6 


— 27:49:22.61 


1.900 


28.25 


1600 


60112 


8.0 ± 1.6 


MHCUloD — 252 


01:58:33.6 


— 24:59:31.10 


2.016 


28.46 


1600 


60112 


5.1 ± 1.3 


FKsll38 — 262 


11:40:48.6 


o/^ o/^ oo r'o 

— 26:29:08.50 


2.156 


29.07 


3000 


17 


8.3 ± 1.6 


TV T T~> "1 n O /I o / ■ o 

MHC1324— 262 


13:26:50.82 


— 26:30:51.10 


2.280 


28.45 


1600 


60112 


1.9 ± 0.8 


Ubbl425— 148 


1/1 r'o r'ri 

14:28:52.50 


"1 01 r'o 

-15:01:37.50 


2.349 


28.66 


1600 


60112 


5.4 ± 1.3 


ijoUoOo VV UUz 


1 V- 1 /I ■ 1 /I 7 


ou.io.zy. iU 


O QOQ 

z.oyo 


Z / . / O 


ooUU 




/I 1 ^ 1 1 

4.1 IIZ 1. i 


MRC0406-244 


04:08:51.5 


-24:18:16.39 


2.427 


29.03 


1600 


60112 


3.8 ± 1.1 


MG2308+0336 


23:08:28.05 


03:36:20.70 


2.457 


28.51 


1600 


60112 


4.5 ± 1.2 


MRC2104-242 


21:06:58.1 


-24:05:11.00 


2.491 


28.84 


1600 


60112 


4.1 ± 1.1 


MRC2139-292 


21:42:16.7 


-28:58:40.00 


2.550 


28.73 


1600 


60112 


4.8 ± 1.2 


MRC0316-257 


03:18:12.0 


-25:35:11.00 


3.130 


28.95 


46000 


3482 


3.8± 1.1 


B20902+34 


09:05:30.1 


34:07:56.89 


3.395 


28.78 


1200 


64 


3.2 ± 1.0 


4C41.17 


06:50:52.1 


41:30:31.00 


3.792 


29.18 


5000 


79 


6.7± 1.5 


TNJ1338-1942 


13:38:26.0 


-19:42:31.00 


4.110 


28.71 


5000 


17 


3.8± 1.1 


6C0140+326 


01:43:43.8 


32:53:49.31 


4.413 


28.73 


5000 


79 


4.5 ± 1.2 



Our analysis focusses on fields around HzRGs with 1.2 < 2 < 3. The other fields are used as control fields. 
^ Radio luminosities from lDe Breuck et al.l 1 12010 ). 



